Quantum field theory as a bilocal statistical field theory 
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We develop a reformulation of the functional integral for bosons in terms of bilocal fields. Cor- 
relation functions correspond to quantum probabilities instead of probability amplitudes. Discrete 
and continuous global symmetries can be treated similar to the usual formalism. Situations where 
the formalism can be interpreted in terms of a statistical field theory in Minkowski space are charac- 
terized by violations of unitarity at very large momentum scales. Renormalization group equations 
suggest that unitarity can be essentially restored by strong fiuctuation effects. 
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I. INTRODUCTION 



The question whether the fundamental constituents of 
nature are particles or waves is discussed since a long 
time. The central equation of nonrelativistic quantum 
mechanics - Schrodingers equation [l[ - is an equation 
for a wave function il)(x,t). On the other side, as was 
first pointed out by Born [2| , this wave function, or more 
precisely its modulus square \'4!{x, can be interpreted 
as the probability density to detect a particle at the point 
X at time t. Indeed, the signals found with experimen- 
tal devises such as fog chambers or silicon detectors are 
localized tracks as expected for pointlike particles. In 
this sense an electron for example has properties of both 
a particle and a field. An alternative view is that they 
always behave like waves and that the particle-like prop- 
erties arise from the special case of a wavefunction that is 
localized in position space, i. e. ip{x,t) ~ S^^^x — x'{t)). 
However, such localized states are usually destroyed by 
the unitary time evolution corresponding to Schrodingers 
equation. To explain the experimentally found particle 
track one has to assume a frequent "collapse of the wave- 
function" occurring in the measurement process. The 
time evolution during the measurement is therefore not 
described by a unitary evolution operator, in contrast to 
situations where no measurement is performed. More- 
over, the fact that this collapse of the wavefunction hap- 
pens instantaneously has lead to interpretational prob- 
lems expressed in terms of different "paradoxa" , most 
prominently the one of Einstein, Podolski and Rosen Q . 

In our modern understanding quantum mechanics is 
the non-relativistic limit of the underlying quantum the- 
ory of fields. In this framework particles are seen as "exci- 
tations of the vacuum" . However, this notion is not very 
precise and seems to be motivated mainly by perturba- 
tion theory for vanishing or small interactions where cre- 
ation and annihilation operators for particles are intro- 
duced similar to the creation and annihilation operators 
for excitations of the harmonic oscillator in quantum me- 
chanics. Concerning the interpretational difficulties men- 
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tioned above, quantum field theory does not lead to much 
new insight. The formulation of the theory is mainly in 
terms of wave functions or fields, but one still needs the 
interpretation in terms of probabilities to "find a particle 
in a certain state" to bring the formalism into contact 
with experiments. This includes also the collapse of the 
wavefunction connected with all the conceptual difficul- 
ties. For a recent, clear written discussion of these issues 
see 3. 

Despite some problems with interpretation or founda- 
tional issues, quantum field theory is a very successful 
theory in praxis. It describes a wide range of physical 
phenomena ranging from the Standard model of elemen- 
tary particle physics to effective theories at low ener- 
gies and momenta. The quantitative agreement between 
some theoretical predictions and experimental observa- 
tion is amazing. It is therefore sensible to investigate 
foundational questions concerning the nature of particles 
(or fields) in this framework. 

A modern formulation of quantum field theory that 
is not restricted to perturbation theory is the functional 
integral approach [5| . This formalism is quite close to for- 
mulations of statistical field theory ; ■ Tlie far reach- 
ing parallels in both formalisms were important for many 
theoretical achievements, for example the development 
of the renormalization group theory 0, One might 
ask whether this duality between quantum and statisti- 
cal field theory has a deeper physical origin and whether 
it can help us for a better understanding of foundational 
questions and of what "particles" are. One might hope 
that some difficult points in the interpretation of quan- 
tum mechanics related to the "collapse of the wave func- 
tion" could be clarified. 

On first sight the functional integral formulation of 
quantum field theory has the disadvantage, that the con- 
cept of particles does not appear explicitly. In contrast 
to the canonical or operator formalism, no creation and 
annihilation operators appear. In fact, the formulation is 
in terms of a weighted sum over field configurations, ex- 
pectation values, correlation functions, fluctuations and 
so on. It is a priori not clear how particle-like prop- 
erties are embedded in this framework. To relate the 
formalism to experiments one has to go an indirect way. 
The results calculated in the functional integral formal- 
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ism must be interpreted by translating them to the op- 
erator language. In this way one can for example extract 
information about the scattering matrix that appears in 
the canonical formalism from the correlation functions 
calculated with the functional integral. Although this 
works in praxis, it would be desirable to have a closed 
interpretation of the functional integral for conceptual 
reasons. This approach has many advantages, for exam- 
ple it is manifestly Lorentz-invariant, can be applied also 
for gauge theories and strong interactions and captures 
non-perturbative features. 

The fact that particles are somehow hidden in the func- 
tional integral formulation can also be seen as a chance. 
After all, also the connection of the canonical formalism 
to experiments needs nontrivial interpretational steps. 
One has to interpret the modulus square of some wave- 
function as a probability density to "detect a particle" 
and the formalism includes the collapse of the wavcfunc- 
tion with all its difficulties. One might hope, that a direct 
interpretation of the functional integral - if it can be es- 
tablished - does not suffer from these problems. One rea- 
son is that the formalism is completely Lorentz-invariant. 
This also being true on the interpretational level could 
solve some problems coming from the fact that the col- 
lapse of the wavefunction by being instantaneous breaks 
Lorentz-invariance. 

In the canonical formalism the physics of n (conserved) 
particles is described by a n-particle state in Fock space 
and its time evolution. Fock states do not have a di- 
rect correspondence in the functional integral formalism. 
Here, the same physics is described by a 2n-point cor- 
relation function. Instead of describing the state at a 
given time, these correlation functions address directly 
the evolution and have two space-time arguments for ev- 
ery particle - typically one corresponding to the initial 
and the other to the final time. At the first rough glance 
this seems to be already close to the experimental situ- 
ation where indeed some sort of correlation between ini- 
tial "measurements" (or preparation of state) and final 
measurements are considered. However, one must keep 
in mind that the correlation functions used in the func- 
tional integral description are not quite the ones known 
from statistical theories. Both meaning and properties 
differ. 

In statistical field theory one calculates expectation 
values, correlation functions etc. by taking the sum over 
all possible field configurations tf weighted by the real 
and positive semi-definite probability measure 

pM=e-^M, (1) 

which we call "functional probability" . For a local theory 
in d dimensions the action S[lp\ is given as an integral over 
the Lagrange density 

Sy]^ j df^xC (2) 

where £ is a function of Lp and its derivatives. In quan- 
tum field theory, the corresponding measure under the 



functional integral 

is complex and can therefore not be interpreted as a prob- 
ability. In addition, to calculate transition probabilities 
from the correlation functions one has to take the modu- 
lus square at some point. Despite this subtleties, one can 
transfer a lot of knowledge and intuition from statistical 
to quantum field theories. 

We propagate here a scenario where the fundamental 
constituents of nature are bifocal microscopic quantum 
fields. In the language of Bell, they are the "beables" 
[lOj . i. e. they exist in the sense of an objective real- 
ity. However, this does not mean, that this microscopic 
fields can be measured directly. Only parts are accessi- 
ble, namely the ones that have particle-like properties. 
Particles are not fundamental objects but arise on a sec- 
ondary, effective or derived level. It is conceivable that 
a description in particle-related terms makes sense un- 
der special circumstances, only. For example, it may be 
that the concept of particles is useful in measurement-like 
situations while in other situations this is not the case. 

The formulation of the field theory is in terms of prob- 
abilities for microscopic field configurations. No assump- 
tion is made whether these probabilities arise from a un- 
derlying deterministic theory. It may well be that this is 
not the case such that the description of nature is prob- 
abilistic already on the fundamental level. 

The idea to abandon particles as the fundamental con- 
stituents of nature was also put forward by Davies . 
From an investigation of the particle concept in curved 
space-time he comes to the conclusion "The concept of 
particles is purely an idealized model of some utility in 
flat space quantum field theory. Away from that limited 
context, however, the concept becomes much less useful 
[...]". Many problems discussed in ref. [ll| might find a 
solution in the scenario proposed here, where the funda- 
mental description is in terms of the microscopic fields 
and their correlations. 

Since in our scenario particles are not fundamental but 
emerging objects, one cannot assign them properties in 
the sense of an objective reality. Partly due to this fea- 
ture, our scenario is in accordance with the theorems of 
Bell [13] or Kochen and Specker Anyhow, the ap- 

proach put forward in this paper does not constitute a 
hidden variable theory, since the description of the funda- 
mental entities - the microscopic fields - is probabilistic. 

This paper is organized as follows. In section [H] we re- 
view briefiy the functional integral formulation of quan- 
tum field theory for bosonic fields and discuss a rewriting 
in terms of a bilocal field. Discrete and continuous sym- 
metries of this approach are discussed in section IIIII and 
positive correlations in section|TVl In the subsequent sec- 
tion|V]we propose a scenario that allows an interpretation 
of the bilocal functional integral formalism in terms of a 
statistical field theory for Minkowski signature. We pro- 
vide arguments for this by investigating renormalization 
group equations for an Euclidean theory with complex 
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action. In section I VII we discuss the properties of a bilo- 
cal statistical field theory and sketch how it can be used 
to describe concrete experiments. Finally, we draw some 
conclusions in section IVTIl 



II. A BILOCAL FUNCTIONAL INTEGRAL 

We start this section with a brief review of the stan- 
dard functional integral for a bosonic scalar field ip(t, x). 
For simplicity, we assume that if has a vanishing vacuum 
expectation value for vanishing external source J. If this 
was not the case, we could easily redefine the fields ac- 
cording to ip'{x) = ip{x) — ipvEV- A crucial object for the 
functional integral formulation is the partition function 

Z[J] ^ I Dip e^^^'^^"'' ^^''''''^^^'^'^^^^'^''^^^'"^^^^ . (4) 



For a local quantum field theory, the action S is given by 
SM ^ [ C, (5) 



where £ is a Lagrange density that depends on the field 
if and its derivatives at the point x. The space-time coor- 
dinate X — {xq , x) contains both the time xq and the spa- 
tial position X and is an abbreviation for J cLxq J (Px. 
Form the partition function we can obtain the source- 
dependent expectation value according to 



{^(x)) 



I 



Z[J] 
1 
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' Z[J] SJ*{x 



■Z[J], 



(6) 



and the assumption made above implies that this van- 
ishes for J = 0. It is useful to introduce the functional 
W[J] according to 



such that 



Z[J] 



{^{x)) ^ 



-iW[.J] 



SJ*{x 



■W[J] 



(7) 



(8) 



Another important object is the two-point correlation 
function for J = 



DFix,y) = {ip{x)ip* (y)) 
which is obtained from the partition function as 
o 1 5 S 



(9) 



DF{x,y) 



Z[Q\ 5J*{x) 6J{y) 
6 S I 
''SJ*{x) SJ{y)^^'^^^J=" 



(10) 



In the last line we used that SW/SJ vanishes for J = 0. 



The two-point function can be related directly to a 
quantity that is accessible to experiments. For a single 
particle that has been detected at the space-time point y, 
the probability to find it at the position x at time xq > ya 
is proportional to 

p{x,y) = \D{x,y)\^ = (11) 

This transition probability can be expressed as a "dou- 
ble" functional integral 



Pix, y) 



1 



Z[0]Z*[0] 



Dip+Dip^ ^P-^-{x)lp-{x) 



x(^+(y)^_(y))*e*^[^+le-^^'['^-l. (12) 

Here we work with the fields ip^ and p- which should be 
considered as independent. Indeed, we can easily pull the 
two functional integrals in Eq. (fT2)) apart, such that the 
terms involving p^ give Dp{x, y) while the terms involv- 
ing ip- give the complex conjugate expression. For usual 
quantum field theories the action S is real up to a small 
term that ensures the convergence of the functional inte- 
gral and is important at a later stage for the description 
of contour integrals. For example, the action S might be 
of the form (we use the metric 77 = diag(— 1, 1,1,1)) 

S[p] = I (p*(x)(a^a^ - m2 + ie)p{x) + 5i„t. (13) 

J X 

The conjugate action S*[(p\ is then of the same form 
with ie replaced by —ie. More generally it is defined 
hy {S[p>]Y ^ S*[p>*]. 

We can obtain the transition probability in Eq. 
also from the generalized partition function 

Z[J+,JJ\ = j Dp>+Dp>^ eyiY>[iS[ip+\-iS*[ipJ\ 



i / {J^p+ + c.c.} + i / {Jlp- + c.c.} 



(14) 

We note that Eq. ([14]) contains the functional integral 
over the field configurations ip^ and One may also 
write this as a single functional integral over the fields 
that depend in addition to the position variable x on the 
contour time tc which is integrated along the Keldysh 
contour (T^ . However, for our purpose it will be more 
convenient to work directly with the expression in Eq. 

It is clear that the modulus square of correlation func- 
tions such as the transition probability p{x, y) can be ob- 
tained from the partition function in Eq. (jl4p by taking 
appropriate derivatives with respect to the sources J+ 
and J_. Since we are interested in the modulus square, 
derivatives with respect to and J_ will allways appear 
in pairs. For example 



p{x,y) 



Z[J+,J-] 
_1 S S 6 

5Jl{x)5J*_{x)5J+{y)5J-{y) ^ +' -^\j+=J-=o 



(15) 
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Such correlation functions can also be obtained from 
a partition function involving bilocal source terms. It is 
defined as 



with 



Z[K] 



J Dipj^Dtf^ explyiS[Lp+\ — iS*[LpJ\ 



{K*{x,y)^+{x)^-{y) + c.c.}\. (16) 



The source field K depends on the space-time variables 
X and y. The transition probability is obtained from 
functional differentiation according to 



p{x,y) 



1 



Z[K] 6K*{x,x) 6K{y,y) 



Z[K] 



K=0- 



(17) 



At this point two remarks are in order. 

(1) The bilocal source term K{x,y) is of bosonic na- 
ture. For bosonic fields ip+ and if- this is obvious but 
also when ip^ and ip- are fermionic fields, it is the case, 
indeed, the source K* couples to the bilinear function 
(p+ip- which is a bosonic operator also for Grassmann- 
valued fields (p+ and 93 „. 

(2) Derivatives of Z[K] at K = always factorize into 
a functional integral over the field (p+ and over the field 
ip^. This is not the case for K ^ 0, however. 

We can formally rewrite Eq. ([T6l) in terms of a func- 
tional integral over a bilocal field xi^TU)- For this pur- 
pose we introduce under the functional integral on the 
right hand side of Eq. a factor of 1 in the form 

1- / DxS[x], (18) 

where S[x\ is a functional Dirac distribution defined by 

Dx S[x]F[x] = F[0] (19) 

for a regular functional of the bilocal field xi^^ u) s-iid its 
complex conjugate. The notation is such that 5[x\ is a 
distribution for both the real and imaginary part of x- 
An example for a concrete realization is 

5[x] = lim c(a)exp ( -a / x*{x,y)x{x,y) ] , (20) 

V Jx,y J 

where c{a) is an irrelevant constant. 

After a shift in the integration variable Xj obtain 
from Eq. (HH) 



Z[K] = J D^+D^^Dx S[x - ip+ip- 



xe 



iS[v+]-iS'[v-]^f^_JK'{x,y)xix,v)+c.c.} 



(21) 



We can now formally perform the functional integral over 
ip^ and ip- such that this becomes 



Z[K] = / Dx 



;-Q[x] ^I^,y{K''{x,v)x{x,y)+x'{x,y)K{x,y)} 

(22) 



Eq. (|22p gives a description of our quantum field theory 
in terms of a bilocal bosonic field x- 



III. SYMMETRIES 

In the following we discuss the transformation behavior 
of Q[x] under various discrete and continuous transfor- 
mations of the argument x- As before we assume that 
(/5+ and ip- are bosonic scalar fields. 

For the bilocal field x(a^, y) the charge conjugate field 
is defined by Xc{x,y) — x*{x,y). For the bilocal action 
Q[x\ we find 

e-Q[xd = I D^+D^^5[x-{ev+){iV-)] (24) 

(25) 



^tS[v+]^-iS-'[v-]^ 



Here ^ is a complex phase (^*^ = 1) that determines 
the intrinsic charge parity of the scalar field After 
relabeling the integration variable ^^p+ — ^ </9+, ~^ 
(/5+, C*<(5- — 'p*_, C^l^ ~^ y^- this gives 



J Dip+Dip^ S[x 



(p+ip-\e 



(26) 

with Sc[^] — S[^*ip*] the charge-conjugated microscopic 
action. It follows that Q[xc] = Q[x] when S is invariant 
under charge conjugation, i.e. S[ip] = Self]- 

Very similar to charge conjugation are parity and time 
reversal. The transformations are defined for the bilo- 
cal field X similar as for the local fields ip+,ip- but as 
for charge conjugation no complex phase (intrinsic par- 
ity) appears for x- In general, Q will show the same 
symmetries as does S, in particular under the combined 
transformation CPT. 

As a final discrete transformation we discuss hermitian 
conjugation of x, i-e. complex conjugation together with 
an interchange of the arguments 



xA^^y) = x*{y,x). 

This gives for the bilocal action 



(27) 



-Qlxs 



j Dip+Dip^ 6[x - (p*^'P*+] e^^['^+le-''^*['^-l. 

(28) 



After relabeling the integration variables according to 

95+ ^ V?!, V?+ t/?-, ip- 95+, V'- ~^ 'P+ this be- 
comes 



nxs 



iS'lvl]piS[ipt] 



Dip+Dip_ S[x - ip+^-] ^^+'e 



(29) 



5 



Note that this imphes that e"'^!-^! is real when restricted 
to hermitian bilocal fields with x — Xs- 

We now turn to the continuous symmetries. As an 
example we consider a global U{1) transformation 

X(x,2/)^e*"x(a;,y)- (30) 
The bilocal action transforms according to 



such that 



= J Dip+Dip^ 6[x- e '"ip+ip- 



g«'S'[¥'+]g-iS'*[v-]_ 



(31) 



For an action S[ip] that is invariant under a global U{1) 
transformation we can change the integration variable, 
for example according to 



which gives 



(32) 



(33) 



In a similar way, other continuous symmetries of S find 
their corresponding counterpart for the bilocal action Q. 



IV. POSITIVE CORRELATIONS 

The fact that Q can be written in terms of a double 
functional integral as in Eq. (j23p has some interesting 
consequences. Although the weight e"*^!-^' in Eq. ((22)) is 
in general complex, some expectation values and correla- 
tion functions are real and positive. Consider for example 
the transition probability in Eq. (jl7p . 



Z[K = 0] 



Dxx{x,x)x*{y.y)e-'^^^l (34) 



Since this can be written as a squared amplitude corre- 
lation function as in Eq. (|12p it is a real and positive 
semidefinite function of the arguments x and y. From 
the expression in Eq. (IMl) this is not obvious but must 
be seen as a property of the bilocal action Q[x\- 

Similar to the transition probability in position space 
one can construct other positive correlation functions. 
Let us consider a complete orthonormal set of functions 
fn(x) with the properties 



(35) 



where 5(x — y) = 5{xa — yQ)5^^\x — y). We define the 
bilocal fields in this basis by 

X{n,m)= f^{x)frn{y)x{x,y), (36) 

J x.y 



x{x,y) = Vx('^,™)/n(a;)/*^(y)• 



(37) 



Note that x(^: is hermitian (antihermitian) with re- 
spect to the arguments n and m if and only if xi^j v) is 
hermitian (antihermitian). 

Correlation functions that involve the modes x("j '^) 
with equal indices n = m will be real and positive semi- 
definite, for example 



Z[K = 0] 



Dxx{n, n)x* (m, ??i)e~^ W > 0. (38) 



This is immediately clear since these functions can be 
written as a product of an amplitude correlation function 
with its complex conjugate. Moreover, we can choose a 
different basis for the field which is labeled with the index 



xin, in) 



9l{x)gfn{y)x{x,y) 



(39) 



(and similar for x*(n, m)) with an orthonormal set of 
functions gfi. We see that the modes of x can be 
written as x("-i ^) in some basis have the special and 
interesting property that expectation values and correla- 
tion functions are positive semi-definite. In the follow- 
ing, these modes are therefore called "positive correlated 
modes" . Of special interest for experiments are such sets 
of functions that are "local in time" , i.e. of the form 



f{t,n){x) =S{xq~ t)(j)n(x) 



(40) 



We combined here the time t and the index n to a com- 
bined index {t,n). The bilocal field in this basis reads 
x((ii, n-i), {t2, n.2)) and expectation values such as 

I J ^XX((i2,n2),(i2,n2))x*((ii,ni),(ii,^ii))e-'3W 

(41) 

describe the probability for a particle that was prepared 
in the state labeled by ni at time ti to be found in the 
state corresponding to the index n2 at time ^2 ■ From the 
equivalence of this to a squared amplitude correlation 
function it is clear that this is a real and positive semi- 
definite function. 

We note that a basis of the form in Eq. pO|) has an 
interesting feature in situations with translational invari- 
ance in time direction of the action S[(p]. In that case 
correlation functions as in Eq. (j4ip are invariant with re- 
spect to an overall shift in time but also with respect to 
a transformation where only the first time argument of 
every bilocal field is shifted, e.g. 



X((t,n),(i,n)) ^X((i + A,n),(i,n)). 



(42) 



This is also clear from the prescription in terms of a 
squared amplitude correlation function. 
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To close this section we remark that the positivity 
property as in Eq. ([551) is independent of the microscopic 
action S[(p] that goes into the definition of the bilocal 
action in Eq. (|23p . It may be a real functional but can 
contain also an imaginary part of arbitrary size. 



V. PROBABILISTIC INTERPRETATION AND 
RENORMALIZATION GROUP FLOW 

So far, the formalism developed in the previous sec- 
tions gives only a reformulation of ordinary quantum field 
theory. Similar to the action S[ip], there is no probabilis- 
tic interpretation of Q[x] since it is in general a complex 
functional and e"*^'-^] cannot be seen as a functional prob- 
ability. (For hermitian x one has e-'3W G M.) We wih 
now examine the meaning of Q[x] more closely and pro- 
pose a scenario that leads to real and positive functional 
probabilities e~'5[^l. 

The bilocal action Q [x] is hard to calculate. Its defini- 
tion in Eq. ()23|) shows that it is in a non-trivial relation 
to the action S[ip] even when the latter is of a Gaus- 
sian form. This makes it hard to develop an intuition 
or to calculate Q[x\ even approximately. The situation 
is better for the partition function Z[K] or the correla- 
tion functions that follow as functional derivatives of this 
object. The correlation functions that can be obtained 
from Z[K] by taking derivatives with respect to K and 
setting K = subsequently can also be obtained as a 
product of two amplitude correlation functions. These in 
turn can be calculated exactly for Gaussian actions S[ip\ 
and many approximation methods exist for theories with 
interactions. 

Often, the correlation functions are directly related 
to objects that can be measured in experiments such 
as scattering cross-sections or decay rates. While cor- 
relation functions are accessible to experiments, this is 
different for the action S[ip] itself. In principle, the com- 
plete knowledge of all (independent) correlation functions 
would allow the determination of S[ip], but in praxis this 
is not helpful. With low energy experiments one can only 
probe an effective action r[i^] where - loosely speaking 
- high energy fluctuations have been integrated out al- 
ready. The determination of S[(p] itself is not possible. A 
connection between the effective action T[ip\ and the mi- 
croscopic action S[^p] is given by renormalization group 
flow equations. Due to the nonlinear structure of these 
equations, it is in praxis not possible to determine S[ip\ 
uniquely by starting from the low energy action T[^p\ and 
following the flow to the ultraviolet. 

A central question in this context is: Could the micro- 
scopic action S[ip\ be of such a form that the exponen- 
tiated bilocal action e~'^^^^ is real and positive? In that 
case we would have a direct probability interpretation 
of the bilocal functional integral in Eq. (1^^ and quan- 
tum field theory could be seen as a bilocal statistical field 
theory. 

One possibility for this scenario is an imaginary micro- 



scopic action 

S[ip] = z^M, (43) 

with a real function S[(p]. From Eq. it follows that 
Q-Q[x\ would be positive semidefinite since it is a func- 
tional integral with a real and positive weight [l^. Is it 
possible that the microscopic action S[(p] is of the form 
in Eq. without changing the effective action r[(p] 
that can be measured experimentally? We will discuss 
this question in the remainder of this section. We note 
however, that the form in Eq. (|43p is not an necessary 
condition for a real and positive functional probability 
Q-Qlx]^ Por example, S[ip] can always contain a constant 
((/^-independent) real part which would change the overall 
normalization of e~*^[^l , only. 

Although we are ultimately interested in quantum field 
theories with Minkowski signature of space-time, it is of- 
ten useful to study theories with Euclidean signature. 
The reason is that both types of theories can be con- 
nected by analytic continuation. As an example we con- 
sider a theory for a bosonic scalar field ip with action as 
in Eq. (|13p . In momentum space the quadratic part reads 

/ ^* (P) (Po -f-m'^+ ie) if{p) . (44) 
Jp 

The ie term is a reminder that we should use the Feyn- 
man prescription to evaluate loop integrals that contain 
the propagator corresponding to the above action. The 
poles for frequency integration are shifted slightly away 
from the real axis 

Po = ±(VF^ + m2 _ (45) 

Equivalently one can drop the ie term and make a rota- 
tion of the frequency integration contour according to 

Po = e^"P4 (46) 

where p4 is now integrated from — oo to oo. The time is 
replaced by the variable r according to 

t = e-*"T. (47) 

The choice a — > 0+ corresponds to Minkowski signature 
and a ^ 7r/2 to Euclidean signature. For the later choice 
the action in Eq. ([T^ becomes 

S^iSe = I I 'p*{xE)i-dl-W^ + m^)^{xE) 

Jxe 

+iSEAnt, (48) 

with xe = (t, x) . The analytic continuation has two 
important features. First, the inverse microscopic prop- 
agator is now of a form where frequency and spatial mo- 
mentum appear in a unified way and with equal sign. 
Second, the functional integral with oscillating weight 
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becomes now a functional integral with the real and pos- 
itive weight 



-Se[v] 



(50) 



The Euclidean action Se is obtained by replacing t by 
—ir at various places. It conventional to add one addi- 
tional overall sign, as well. 

Note that a functional integral with real and positive 
weight e*"^!"^! = e^'^^'^l in Minkowski space becomes after 
analytic continuation to the Euclidean space a functional 
integral over the complex and oscillating weight 



-Se[v>] 



-iSeIv] 



(51) 



To summarize, an complex oscillating weight functional 
with Minkowski signature corresponds to a real and pos- 
itive weight functional with Euclidean signature, while a 
real and positive weight functional with Minkowski sig- 
nature corresponds to a complex oscillating weight func- 
tional with Euclidean signature. 

Of course, it is known from many experiments that the 
actual world has Minkowski signature. On the other side, 
the investigation of field theories in a functional integral 
is simpler with Euclidean signature. For example, the 
construction of regulator functions is straightforward in 
the later case and exact renormalization group equations 
can be used directly. 

Let us consider the flowing action rfe[(^]. It is defined 
as the effective action in presence of an infrared cutoff 
term 



(52) 



in the microscopic action. More precisely, we modify the 
Euclidean version of the functional integral in Eq. ^ 
according to 



(53) 

with AS**; as in Eq. ([52]) . Consider now the Legendre 
transform of Wk [J] 



ffe[(^] = / {J*ip + c.c.}-Wk[J], 



(54) 



with tf ~ SWk[J]/SJ. The effective average action or 
flowing action is defined by subtracting from this the cut- 
off term 



Tk[ip] = fk[ip] - ASk 



(55) 



The flowing action interpolates between the micro- 
scopic action 



lim Tk[ip] = S[(p] 

k^oo 



(56) 



for large values of the infrared cutoff parameter k and 
the quantum effective action 



lim rfe[v3] = r[ip] 

fe— >o 



(57) 



for vanishing infrared parameter k. We used here the 
properties of the cutoff function Rk{q^) — ?> for fc — > 



and Rkiq^) 



for fc — i- CX3. The quantum effective 



action is the generating functional of one-particle irre- 
ducible Feynman diagrams. For nonzero value of fc, the 
flowing action FfcM is closely connected to an action for 
averages of fields flSI]- 

With these definitions, the flowin g ac tion obeys the 



exact functional differential equation 16 1 



dkTk 



(58) 



The trace operation Tr includes an integral over momenta 
and a sum over possible internal degrees of freedom. We 
note that the right hand side of Eq. (|58p contains F^ ' , 
the second functional derivative of F^ [tp] 



6ip* {q)6(f{q') ' 



(59) 



Since both sides have a functional dependence on ip, Eq. 
()58p is a functional differential equation. For an review 
of the flow equation and the concept of the flowing action 
we refer to ref. [l7j 

We will now investigate the renormalization flow of an 
action that is rotated from the real axis to the complex 
plane. In the microscopic regime for large k it is of the 
form 



(60) 



with Ffe[(/3] a real Euclidean functional of the held if. 
We choose the parameterization such that afc = cor- 
responds to the usual sign convention for an action with 
Euclidean signature. It is convenient to choose the cutoff 
function according to 

ASk = e''"' f V*{p)Rk{p'Mp) (61) 
Jp 

with real and positive Rk{p^)- The flow equation in Eq. 
becomes now 

e"'''{dkTk+i{dkak)Tk) 

= ^Tr(fi'^ + Rk)-\dkRk + i{dkak)Rk). (62) 

This should be seen as a complex equation that deter- 
mines both SfeFfe and dkUk- 

In general, the simple form in Eq. (j60p will not be 
conserved by the flow. In an expansion of the flowing 
action 



(63) 



with real operators ©'"-"[(/j], the coefficients 7^"'' will be 
complex functions of the scale parameter k. In princi- 
ple one can write 7^"-' = e'"*! t-^""* with 7^"'' £ R and 
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investigate the flow equations for the set of phases a^" . 
For simphcity we assume here that all these phases are 
equal, a^"^ = a/t, so that we can use Eq. (15^ . This 
gives useful qualitative insights. For concrete models one 
can of course also study truncations with a number of 
independent complex coefficients 7^"''. 

We note that both sides of equation (|62|) are real for 
Q!fc = such that this corresponds to a renormalization 
group fixed point. Moreover, there can be only fixed 
points for ak when e'"'' is real, since the right hand side 
of Eq. (|62|) is real for dtak — 0. In particular, the flow al- 
ways drives the action away from the point with ak — vr/2 
corresponding to Eq. ()51|) . In situations with strong in- 
teraction effects, the renormalization group modiflcations 
are large. In this case, the action F^ [</)] or more concrete 
the expansion coefficients 7^"'' in Eq. ()63|) will be driven 
towards the real axis. 

These arguments show that ak = is likely to be an 
infrared attractive renormalization group flxed point for 
a large class of initial microscopic actions S[ip]. In sit- 
uations with strong running due to large fluctuation ef- 
fects, this fixed point is approached very quickly. By this 
mechanism, renormalization group modifications could 
render a theory with complex or purely imaginary ac- 
tion S'£;['p] = iSsl^] Euclidean signature to a real ac- 
tion TeI^] on a macroscopic scale. Translated back to 
Minkowski signature, this implies that for strong inter- 
actions, an imaginary action S[(j)] = iS[ip\ corresponding 
to a positive weight 

could lead to an effective action T[(p] of the usual form. 

This opens the possibility for a probability interpreta- 
tion of quantum field theory in terms of the bilocal action 
Q[x]- Before we discuss this in more detail in the follow- 
ing section, let us make some remarks on the connection 
of the formalism to usual quantum field theory. 

(1) For free theories where F^ and higher terms van- 

ish, the inverse propagator F^ is not renormalized. This 
implies in particular that its complex phase remains con- 
stant so that the above mechanism is not applicable. In- 
teraction effects are crucial in this respect. 

(2) One might wonder, whether the proposed scenario 
can be applied to the Standard model, where most inter- 
action parameters are rather weak. However, although 
one needs a strong renormalization group running, it is 
nevertheless possible that the infrared physics is domi- 
nated by weak coupling constants. The reason is that in 
most cases strong interactions are strongly screened by 
quantum fluctuations leading to effective theories with 
small interaction parameters in the infrared. It may 
be that the strong renormalization group running which 

drives the coeflacients 7^."' towards the real axis happens 
at very large scales, possibly above all other relevant 
scales such as the electroweak scale or the Planck scale. 

(3) We note that a given (real) effective action F[0] 
on the macroscopic scale can follow from several micro- 



scopic actions. In particular, since the theory-subspace 
with ak = corresponds to a fixed point for ak, there 
will usually be a trajectory in this subspace leading to a 
particular real effective action F[(/)] at fc = 0. 

(4) Usually one assumes that the effective action F[(^] 
in Minkowski space has to be real since this corresponds 
to unitarity of the S-matrix. Euclidean quantum field 
theories must have Osterwalder-Schrader reflection posi- 
tivity which implies positivity of the norm for Minkowski 
signature [ll]. It is clear that the scenario proposed 
here leads to deviations from positivity in the micro- 
scopic regime. To prevent conflicts with well established 
and experimentally tested features of quantum mechan- 
ics, unitarity or Osterwalder-Schrader positivity must be 
effectively restored at macroscopic scales. More detailed 
investigations for concrete models must reveal that this 
is actually the case. 

(5) Finally, we remark that deviations from unitarity in 
the microscopic regime could lead to measureable effects 
that could be tested in experiments. Presumably this 
is not possible if unitarity is violated only at very large 
momentum scales but more detailed investigations are 
necessary to clarify this point, as well. 

VI. PROPERTIES OF BILOCAL STATISTICAL 
FIELD THEORY 

In this section we assume for the time being that a 
mechanism as the one proposed in the previous section 
allows to reformulate quantum fleld theory in terms of a 
bilocal statistical fleld theory. We investigate some con- 
sequences and discuss how known physics emerges. 

On the flrst sight, statistical field theory for bilocal 
fields and with a functional probability 

e-QW (65) 

is very different from the quantum mechanical descrip- 
tion we are used to. For example, Eq. (|55|) describes 
probabilities for fleld configurations in space-time, while 
quantum mechanics describes probabilities for measuring 
particles in a certain state conditional to some prepara- 
tion experiment in the past. This is irritating if there 
is only one world or only one space-time so that such a 
probability does not seem to make much sense. How- 
ever, one can interpret Eq. (j65p as describing probabili- 
ties for microscopic field configurations which cannot be 
directly accessed by measurements. We will argue below 
that only particular modes correspond to particles while 
others should be seen as part of a non-trivial vacuum. 
Non-Gaussian terms in Q[x] lead to connections between 
these two sectors. It is natural to employ a statistical 
ensemble for the microscopic or vacuum part of x- 

Most fluctuations of the bilocal field will decay on short 
time scales or the values of x will be uncorrelated after 
some relatively short time. Such fluctuations correspond 
to the vacuum. A typical particle detector produces sig- 
nals only if an excitation is coherent for a long enough 
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time. For a statistical field theory the relevant object is 
the two-point correlation function. For a bilocal theory 
it reads 

G{x+,x-;y+,y-) 

= -^^^—^ J Dxx{x+,x-)x*{y+,y-)- (66) 

The inverse propagator G^^ is defined by 

/ G"^(x+,a;_;z+,z_)G(z+,z„;?/+,y_) 

= S{x+ -y+)S{x_ - y_). (67) 

Typically, it can be written as 

G^'^{x+,x^;y+,y^) = S{x+ - y+)S{x^ - y-)P+Pl, 

(68) 

where P+ and P_ contain derivatives with respect to 
and X- . Modes of x that correspond to long range fluc- 
tuations are characterized by 

P+Plx{x,y)=0. (69) 

Short range fluctuations have a certain influence on the 
form of G~^, in particular they change the values of 
appearing constants such as the mass or coefficients of 
derivative terms. This are the well known renormaliza- 
tion group modifications. 

Equation (|69p cannot be a sufficient condition for a 
mode of x to be directly measurable in experiments. The 
field X still has two independent variables x and y. This 
is in contrast to our experience of a "local world" where 
one can only measure fields depending on a single space- 
time coordinate. 

One can make progress by assuming that a mode that 
is detectable needs not only long-range (or long-time) 
correlations but also positive ones. By this we mean that 
the two-point correlation function G needs to be real and 
positive semi-definite. Modes that have this property are 
the "positive correlated modes" discussed in section ITVl 
To discuss the measurement of particles in a particular 
basis 0„ (x) it is useful to write the field x £^t equal time 
argument as 

Xiit,x),{t,y)) = ^x((i,«),(i,"))0«(x)0:(y) (70) 

n 

+ X{{t,n),{t,m)) 4>n{x)(t)l,{y). 

The first sum on the right hand side describes the con- 
tribution of positive modes while the second sum adds 
off-diagonal terms. The diagonal terms are the ones that 
describe particles according to the rules of quantum me- 
chanics. We emphasize that the fact that the correla- 
tions of x((^j n-)j {ty n-)) are real and positive semi-definite 
is a special property of the functional probabilities e"'^!-^! 
that can be written as a double functional integral as in 



Eq. (|28p . For more general bilocal theories this will not 
be the case. 

We note that this can explain why our world is "to 
large extend local" although the field x is bilocal. The 
"positive correlated modes" x((^j (^i "-)) have only one 
free time argument remaining. The spatial dependence 
is encoded by a function (j)n{x) as in Eq. ()40|) . In many 
cases this function will be localized in some volume of 
space and we encounter the same kind (and degree) of 
nonlocality as in quantum mechanics. 

A useful physical picture is that the fields that are 
relevant for experiments are coarse-grained or averaged 
over some time interval At. The coarse-grained fields 
might be given by 

Xriit,x),it,y)) = Ax{it,S),it,y)) (71) 

where a{x) is a positive function that vanishes for |a:;| ^ 1 
and has the property J a{x)dx = 1. If a mode of x 
is oscillating with frequency large compared to 1/At it 
will be effectively annihilated by the operator A. Os- 
cillating modes typically correspond to negative or more 
general complex two-point correlation function. In con- 
trast, modes where G is real and positive are usually 
constant and therefore non annihilated by the coarse- 
graining. Note however that also the diagonal part of 
X, i.e. the flrst term on the right hand side of Eq. (|70l) . 
might contain terms that are annihilated by the coarse- 
graining in a concrete situation. 

Let us now discuss some physical situations to find out 
whether the picture that emerged so far makes sense. We 
first consider a situation where a particle has been pre- 
pared in an eigenstate of the (single particle) Hamiltonian 
labeled by n = 1 at time t. This situation corresponds 
to the coarse-grained fields 

X((i,l),(t,l)) ^ 0, 

X{{t,n),{t,n)) = fornix. (72) 

The fields at later times t' > t are positive correlated 
to the ones at time t and since n labels eigenstates of 
the Hamiltonian, the cross-correlations vanish. In other 
words we have 

{Xiit',n'),it',n'))x*iit,n),{t,n))) ^ const Sn,n. (73) 

This implies that only the mode with n ~ 1 will be 
nonzero at later times. This situation is precisely as we 
expect it form the laws of quantum mechanics. 

Now let us consider the situation where the measure- 
ment at time t' is performed in another basis which we 
label by the index n. We have now correlation functions 

ixiit' ,fi'),{t' ,fi'))x*{{t,n),it,n))) ^ const pn'n (74) 

that equal by construction the transition probabilities 
from quantum mechanics pa'n (up to a constant factor). 
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Similarly, one can expect that other correlation func- 
tions are equal for quantum field theory and our bilocal 
statistical field theory. However, when one attempts to 
make this more precise, one is confronted with a puz- 
zle: Why does a particle-like excitation not split into 
several ones? Indeed, one could expect nonzero field 
x{{t' , h) , {f , n)) for several n if the correlation function 
in Eq. (|74p with n = 1 is nonzero. The question is closely 
connected with the question why particle number is con- 
served and in particular why it is discrete. For a theory 
with a global U{1) symmetry one can derive conserva- 
tion laws for appropriate currents. However, this does 
not imply the discreteness of the particle number. It may 
well be that this puzzle can be solved once the formal- 
ism is extended to fermionic theories. A certain kind of 
discreteness arises naturally for functional integrals over 
Grassmann numbers. However, this is beyond the scope 
of the present paper and we leave this point open for 
future investigations. 

To close this section we remark that the basis chosen 
to perform the initial preparation at time t does not have 
to be the eigenbasis of the Hamiltonian. It is straightfor- 
ward to transfer the above considerations to more general 
situations. 



VII. CONCLUSIONS 

We have developed a formulation of quantum field the- 
ory for bosonic scalar fields in terms of a functional inte- 
gral over bilocal fields. This has the advantage that phys- 
ical observables such as correlation functions or transition 
probabilities can be obtained directly instead of first cal- 
culating probability amplitudes which need to be squared 
to give probabilities. 

The bilocal formulation can be interpreted in terms 



of a statistical field theory if the bilocal action Q[x] is 
real. An explicit example for this corresponds to a quan- 
tum field theory where unitarity is broken at microscopic 
scales. The flow equation for the effective average ac- 
tion in Euclidean space suggests that is becomes real at 
macroscopic scales even if the microscopic action is imag- 
inary. Strong fluctuation effects can lead to an effective 
action of the usual form at small values of the cutoff pa- 
rameter. In Minkowski space this corresponds to a the- 
ory where the weight of the functional integral is real and 
positive. This implies in particular that the exponenti- 
ated bilocal action 6""^!-^^ is positive. Unitarity would be 
an emerging phenomenon instead of being a property of 
the fundamental theory. 

It is clear that this scenario needs (and deserves) fur- 
ther investigations. A more detailed renormalization 
group analysis of statistical held theories with Minkowski 
signature could reveal interesting relations to correspond- 
ing quantum held theories. By investigating concrete and 
realistic models one could explore the experimental sig- 
natures for deviations from unitarity in the microscopic 
regime 

We have also discussed how properties of particles can 
arise in a bilocal statistical fleld theory. Positive correla- 
tions at subsequent times play an important role. Quan- 
tum mechanical observables can be represented as corre- 
lation functions of bilocal flelds. 

It is planed to extend the formalism to fermionic the- 
ories described by Grassmann valued flelds. The bilocal 
fleld will always be of bosonic nature since it corresponds 
to a pair of original flelds, but one can expect that some 
properties will differ. The formalism in its present form 
already provides evidence for the intriguing possibility 
that quantum fleld theory can be understood as a bilocal 
statistical field theory. 
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